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STUDIES  0?  HnFn  OLIGOMERIC  SPECIES 


Abstract 


Infrared  matirx-isolation  spectra  of  HF  and  HF/H00  mixtures  suspended  in 

'/ 

solid  argon  have  been  obtained  in  the  ranges  400-4000  cm  • Bands  attributed  to 

rl  ' 

HF  and  to  HF-H2O  oligomers  were  observed.  HF  stretching  modes  were  assigned  to 
dimer,  chain  trimer,  cyclic  tetramer,  and  cyclic  hexamer  of  HF.  Some  of  the 
FH*  • *F  angle  deformation  modes  were  observed  for  these  molecules.  The  oligo- 
mers are  formed  during  matrix  deposition  and  the  chain  trimers  appear  to  be 
matrix  stabilized  in  a double  substitution  site  in  argon.  Application  of  a non- 
ideal gas  model  to  the  literature  vapor  phase  data  on  HF  suggests  that  the  dimer 
and  the  hexamer  are  the  dominant  vapor  phase  oligomers,  with  the  hexamer  most 
Important  near  saturation  vapor  pressures.  The  tetramer  is  the  next  most  impor- 
tant vapor  phase  oligomer,  but  it  contributes  only  weakly  to  the  vapor  phase 
properties.  Minimal  basis  set  ST0-3G  computations  on  chain  and  cyclic  HF  trimer 


and  tetramer  molecules  suggest  the  marked  increase  in  stability  that  arises  on 
catenation  or  (unstrained)  cyclization  of  HF  systems 
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STUDIES  OF  HnFn  OLIGOMERIC  SPECIES 
I.  RESEARCH  OBJECTIVES 

It  was  proposed  to  use  matrix-isolation  spectroscopy  to  study  the  properties 
and  vapor  phase  populations  of  oligomers  of  HF  and  DF.  The  properties  to  be  deter- 
mined included  the  existence  of  specific  oligomers,  their  molecular  symmetry  and 
qualitative  geometry,  their  vibrational  frequencies,  and  at  least  a partial  analy- 
sis of  their  molecular  force  fields.  The  research  effort  was  designed  to  provide 
information  concerning  intermolecular  interactions  between  HF  molecules  and,  there- 
by, to  provide  information  that  is  useful  for  modeling  collisions  and  intermolecu- 
lar energy  transfer  in  the  HF  vapor.  Accurate  models  to  describe  collisional 
relaxation  processes  involving  monomer-monomer  collisions  or  monomer-oligomer 
collisions  are  essential  for  providing  quantitative  kinetic  models  for  HF  or  DF 
lasers.  These  lasers  are  useful  for  Air  Force  operations  involving  reconnaissance, 
communication,  detection,  surveillance,  and  weaponry. 

II.  STATUS  OF  RESEARCH  EFFORT 
A.  Matrix-Isolation  Experiments 

Hydrogen  fluoride  has  been  studied  using  matrix-isolation  sampling  several 
times  previously  (1,2,3)  but  the  research  effort  was  directed  at  investigating 
the  behavior  of  matrix-isolated  HF  monomer  rather  than  at  elucidating  properties 
of  H^Fn  oligomers.  The  matrix-isolated  HF  monomer  is  of  interest  because  it 
undergoes  hindered  rotation  in  rare  gas  lattice  trapping  sites.  Rotation  is 
suggested  by  observations  on  several  bands  attributed  to  vibration-rotation  transi- 
tions, and  this  phenomenon  is  consistent  with  results  obtained  on  several  other 
small  molecules  that  are  capable  of  fitting  into  substitutional  trapping  sites  in 

i 

the  rare  ga§  lattices  * e.  g.  H20(4) , NH^(5),  HC1(6),  CH^(7). 


The  only  published  spectra  that  pertain  to  oligomers  are  shown  in 

reference  1,  but  the  HF  samples  shown  there  are  contaminated  with  large  quan- 
tities of  water,  and  many  of  the  observed  absorptions  are  evidently  due  to  anci 

complexes  between  HF  and  1^0  rather  than  to  HnFn  oligomers.  The  spectra  obtained 
in  the  present  research  differ  significantly  from  these  results,  as  water  and 
other  impurity  molecules  are  nearly  all  eliminated.  In  contrast  to  the  status 
concerning  matrix-isolated  HnFn  oligomers,  rather  complete  studies  have  been 
reported  for  hydrogen  chloride  oligomers(8)  and  for  water  oligomers (9) . 

EXPERIMENTAL  TECHNIQUE.  The  matrix-isolation  technique  consists  of  trap- 
ping the  sample  molecules  in  a large  (hundreds-fold  or  thousands-fold)  excess  of 
inert  matrix  material.  In  the  present  research  the  matrix  material  was  exclu- 
sively argon.  Neon  is  probably  a better  matrix  choice;  however,  neon  could  not 
be  used  with  the  available  low  temperature  refrigerator.  The  inert  matrix  inter- 
acts little  with  the  isolated  sample  molecules,  which  therefore  exhibit  near-vapor 
phase  vibrational  band  origins.  The  absorption  peaks  are  sharp  due  to  the 
absence  of  a rotational  band  envelope  or  to  strong  intermolecular  interactions. 

The  rigid  matrix  prevents  coalescence  of  smaller  species  to  form  solid  HF. 

In  "premixed"  experiments  the  gaseous  HF  and  Ar  are  mixed  together  in  a 
sample  cylinder  at  (nominally)  known  pressures.  The  gaseous  mixture  is  then 
sprayed  onto  the  cold  (11  K)  sample  holder,  where  it  freezes  into  a rigid  solid. 

In  the  case  of  HF,  accurate  knowledge  of  the  M/S  (matrix  to  sample)  ratios 
seems  precluded  in  the  work  performed  to  date  because  the  HF  is  strongly  adsorbed 
onto  the  walls  of  the  sample  cylinder  and,  moreover,  it  tends  to  react  with 

these  walls  to  form  SiF,  (stainless  steel  contains  small  amounts  of  Si,  which 

* 

the  HF  manages  to  reach). 


In  this  research  two  types  of  matrix-isolation  sample  preparations  have 
been  performed.  These  are  the  "premixed"  vapor  deposition  just  described  and  a 
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"pin  hole11  or  double  jet  deposition  technique.  In  the  latter , the  HF  vapor 
source  Is  thermos tatted  at  constant  pressure  (the  present  range  of  Interest  is 
less  than  one  atmosphere)  and  then  It  is  allowed  to  escape  into  the  cryostat 
system  through  a tiny  orifice.  The  effusing  HF  beam  is  codeposited  with  an  argon 
beam  on  the  cold  sample  holder.  The  objective  of  such  double  jet  experiments  is 
to  obtain  reliable  spectra  of  the  oligomers  present  in  the  original  HF  vapor  at 
known  temperature  and  pressure.  In  premixed  sample  preparations  the  HF  partial 
pressure  is  too  low  to  allow  the  presence  of  vapor  phase  oligomers  and  they  are 
formed  in  the  matrix- isolation  samples  upon  deposition  by  coalescence  of  monomers 
brought  together — using  concentrated  M/S  ratios.  On  the  other  hand,  hexamers 
and  other  oligomers  can,  in  principle,  be  directly  obtained  from  HF  vapor  at 
very  dilute  M/S  ratios  using  the  double  jet  deposition  system.  Ideally,  their 
absorption  intensities  will  reflect  the  original  HF  vapor  phase  composition. 

In  the  present  deposition  system,  the  pin  hole  is  an  orifice  laser  drilled (10) 
at  the  apex  of  a cone  bored  into  the  outer  end  face  of  a standard  Gyrolok  1/4 
inch  tube  cap.  Orifice  diameters  of  0.005  to  0.0005  inch  were  prepared.  A dozen 
experiments  were  performed  using  double  jet  deposition,  with  results  that  differ 
markedly  from  those  obtained  using  standard  premixed  samples,  but  it  does  not 
appear  that  these  matrix-isolation  samples  reflect  the  vapor  compositions  of 
the  original  vapor.  There  are  still  experimental  problems  associated  with  per- 
fecting the  double  jet  apparatus  but  one  difficulty  may  be  the  pinhole  orifice 
itself.  A more  appropriate  construction  is  to  drill  the  pin  hole  into  a cap  of 
very  thin  foil,  thus  avoiding  the  presence  of  a volume  where  collisional  relaxa- 
tion and  partial  low  pressure  equilibration  can  occur.  The  spectra  to  date 
reflect  a lower  oligomer  composition  than  expected  from  the  original  HF  gas 


, 


» 


sample. 

It  was  originally  intended  to  construct  an  all  monel  vacuum  system  for  the 
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HF  purification  and  sample  preparation.  The  present  system  includes  monel 
tubing,  fittings,  and  valves,  but  also  stainless  steel  valves  and  a 0.5  1 stain- 
less steel  sample  preparation  cylinder.  We  were  forced  to  wait  nine  months 
for  the  sample  preparation  cylinder  and  nearly  as  long  for  some  of  the  other 
components.  The  sample  cylinder  is  the  essential  component  of  the  "premixed" 
matrix- isolation  experiments  and,  in  spite  of  repeated  serious  attempts,  satis- 
factory experiments  were  not  possible  until  the  high  quality  stainless  steel 
vessel  was  finally  delivered.  The  earlier  experiments  used  a Teflon  plastic 
gas  sampling  bag,  which  leaked,  a metal  storage  cylinder  already  present  in  the 
laboratory,  which  reacted  with  HF  to  yield  numerous  unknown  impurities,  and  a 
locally  constructed  stainless  steel  cylinder,  which  reacted  with  the  HF  to  yield 
objectionably  large  amounts  of  SiF^.  At  present  we  have  saved  spectral  charts 
for  above  forty  numbered  matrix-isolation  experiments  (including  some  poor  ones). 
These  experiments  are  for  HF  premixed  with  Ar  at  various  matrix/ sample  ratios 
nominally  ranging  from  about  30  to  several  hundred*  pin  hole  deposits,  and  a series 
of  experiments  for  HF/Ar  with  water  (0-16  and  some  with  0-18)  added  In  order  to 
study  the  spectra  of  the  HF-H^O  complexes  as  well  as  the  HnFn  oligomers.  The 
HF  gas  was  purchased  from  Matheson  Company  and  purified  for  each  experiment  by 
pumping  on  samples  trapped  at  77  K and  then  at  about  185  K.  A partial  spectral 
analysis  has  been  performed  on  the  available  data,  as  discussed  below,  but  before 
publication  additional,  confirmatory  experiments  using  deuterated  molecules,  and 
water  0-18  will  be  performed.  Those  will  allow  force  field  calculations  to  be 
performed  in  some  cases.  The  experiments  are  presently  being  performed  without 
AF0SR  aupport. 

The  research  is  performed  using  an  Air  Products  Displex  Model  202  closed 
cycle  low  temperature  refrigerator  rather  than  the  liquid  helium  cryostat.  This 
is  economical  and  allows  the  use  of  an  argon  matrix,  but  is  not  cold  enough  for 
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a neon  matrix.  A Beckman  IR-9  spectrophotometer  is  used  to  record  the  spectra 
in  the  range  400-4000cm  \ We  hope  to  perform  experiments  in  the  spectral  range 
below  400  cm  ^ using  a Beckman  IR-11  far-infrared  spectrometer  that  is  available 
in  this  laboratory. 

RESULTS.  The  argon  matrix-isolation  spectra  show  intense  absorptions  in 
the  spectral  region  between  3000-4000  cm  ^ that  are  due  to  HF  stretching  modes,  and 
generally  weaker  absorptions  (with  two  exceptions)  in  the  region  400-800  cm  * 
that  are  due  to  bending  or  angle  deformation  modes.  The  region  above  3700  cm  1 
shows  four  intense  absorptions  attributed  to  oligomers,  along  with  a number  of 
weaker  bands.  The  lower  frequency  region  shows  surprisingly  few  absorptions. 
Representative  spectra  are  shown  in  Fig.  1 and  the  absorptions  are  listed, 
along  with  those  of  previous  investigations,  in  Table  I.  These  results  are  dis- 
cussed according  to  the  various  oligomers  in  the  following  paragraphs. 

Monomer.  The  present  experiments  were  not  designed  to  study  the  monomer 
spectrum  which,  as  mentioned,  is  of  interest  in  terms  of  hindered  rotation  of 
HF  in  the  rare  gas  matrix  environment.  The  monomer  absorptions  consist  of  a 
very  intense  doublet  and  two  very  weak  side  bands.  The  weak  side  bands  were  not 
observed  in  the  present  work.  The  intense  doublet  closely  resembles  that  illus- 
trated by  Mason,  Von  Holle  and  Robinson (3),  with  components  of  rather  similar 
intensity  and  not  well  resolved.  The  slightly  more  intense,  higher  frequency 
component  has  been  attributed  to  R(0)  of  the  rotating  monomer;  the  lower  fre- 
quency component  attributed  to  HF  in  a different  matrix  trapping  site. 

Dimer.  The  approximate  geometry  of  the  dimer  is  known  from  microwave  and 
radiofrequency  studies  by  Dyke,  Howard  and  Klemperer (11)  using  molecular  beam 


samples.  The  dimer  has  the  geometry  110-120° 

F • • • H-F 


FF  = 27.9  nm 
FH  = 9.2  nm. 


This  tetratomic  molecule  possesses  six  infrared  active  vibration  frequencies,  but 
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the  only  previous  experimental  observations  consist  of  the  measurement  of  weak 
infrared  absorption  that  occurs  between  several  of  the  intense,  low  J,  vibra- 
tion-rotation bands  of  the  gaseous  monomer.  This  absorption  has  been  studied 
by  Smith(12),  Herget,  Gailar,  Lovell  and  Nielsen(13),  and  Himes  and  Wiggins(14). 
The  absorption  profile  is  poorly  characterized  because  of  experimental  diffi- 
culties, but  there  appear  to  be  three  maxima  and  fine  structure  is  resolved,  over 
portions  of  the  profile.  There  is  an  intense  band  in  the  argon  matrix  at  3881  cm 
that  corresponds  reasonably  well  with  the  center  of  the  complex  vapor  phase  dimer 
absorption.  The  argon  matrix  band  is  assigned  to  dimer  because  it  is  the  last 

oligomer  absorption  to  persist  at  dilute  M/S  ratios,  while  Smith(l2)  assigns  the 

2 

vapor  phase  absorption  to  dimer  because  its  intensity  shows  a P pressure  depen- 
dence. 

The  vibrational  band  origin  may  be  taken  to  be  near  the  388  1 cm  ^ argon 
matrix  value.  The  other  "bands"  observed  in  the  vapor  phase  are  either  hot  bands 
or  else  part  of  a complex  profile  associated  with  the  3881  cm  ^ fundamental.  The 
presence  of  a complex  absorption  profile,  in  contrast  to  that  of  a conventional 
near-symmetric  top  (cf.  Ueda  and  Shimanouchi (15) ) , is  consistent  with  the  avail- 
able knowledge  of  the  dimer.  The  microwave  spectra (11)  show  complexity  suggest- 
ing that  the  molecule  is  nonrigid  and  that  tunneling  of  the  H atoms  occurs: 

H 

^F  • • • H-F  -*■  F-H  • • • F 

Si  ^'H 

Experience  with  other  nonrigid  molecules  has  illustrated  the  complexity  that 
arises  in  their  vapor  phase  infrared  spectra.  For  example,  carbon  suboxide  is 
a quasilinear  molecule  with  a very  low  frequency,  large  amplitude  OOC  angle 
bending  mode  that  is  somewhat  analogous  to  the  rocking  of  the  "free"  HF  moiety 
of  the  H2F2  dimer.  The  appearance  of  the  ^ow  rti80^ut^on  spectrum  is  decep- 


tive(16), as  made  obvious  by  later  high  resolution  studies (17),  and  it  is  only 
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recently  that  the  analysis  of  the  potential  function  and  intramolecular  dynamics 
are  becoming  understood(18) , Likewise,  OH  stretching  modes  of  hydrogen  bonded 
carboxylic  acid  dimers  have  enormously  broad,  complex  band  contours,  even  in  a 
neon  matrix (19) , due  to  intramolecular  interactions  involving  the  hydrogen  bond 
linkages (20) . The  extreme  band  broadness  and  intensification  caused  by  H-bond- 
ing  in  the  carboxylic  acid  dimers  is  not  shown  by  H^F^  dimer  in  the  vapor  or  in 

f 

the  argon  matrix,  but  it  is  present  in  solid  HF  where  large  scale  intermolecular 
couplings  can  occur.  This  is  analogous  to  the  case  for  water. 

The  second  HF  stretching  mode  is  much  less  intense  than  the  3881  cm  1 band, 
and  considerations  of  the  overall  matrix-isolation  data  that  are  presently  avail- 
able suggests  that  the  best  choice  is  a band  at  3787  cm  ^ . The  intensity  of  this 
band  correlates  better  with  that  of  3881  than  do  intensities  of  the  other  possible 
candidates.  Unfortunately,  it  appears  that  matrix  conditions  strongly  affect  the 
intensity  of  the  various  HnFn  oligomer  transitions.  This  is  undoubtedly  because 
of  the  great  flexibility  of  the  oligomer  molecules,  and  it  provides  a problem 
because  a precise  correlation  of  the  band  intensities  with  relative  populations 
is  precluded  as  M/S  ratios  are  varied.  However,  relative  correlations  are  still 
possible  in  many  cases,  as  shown  by  plotting  relative  intensities  against  matrix 
dilutions . 

Assigning  the  3787  cm  ^ band  as  the  second  HF  stretching  mode  of  the  dimer 
is  in  accord  with  published  ab  initio  quantum  mechanical  calculations (21) . These 
place  the  two  HF  stretches  only  43  cm  1 apart  (at  4081  and  4038  cm  *)  although 
they  are  about  200  cm  ^ higher  than  the  observed  values.  In  spite  of  the  hydro- 
gen bond,  the  two  HF  stretches  have  nearly  the  same  vibrational  frequencies,  and 
they  can  be  expected  to  strongly  couple.  Because  of  its  greater  intensity  the 
3881  cm  ^ band  must  correspond  approximately  to  an  "in-phase"  motion  of  the  two 
HF  vibrators,  while  the  weaker  3787  cm  * band  must  correspond  to  "out-of-phase" 
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motion.  The  lower  frequency  mode  has  a larger  contribution  from  the  H-bonded 
HF  moiety  but,  as  noted,  it  does  not  show  the  large  intensif ica tion  and  enor- 
mous frequency  shift  associated  with  some  H-bonded  stretching  modes. 

Two  bending  modes  of  the  dimer  correspond  roughly  to  in-plane  and  out-of- 
plane  motion  of  the  H-bonded  (F*--HF)  H atom.  These  modes  are  calculated  quantum 
mechanically  to  occur  at  588  cm  ^ (in-plane)  and  519  cm  * (out-of-plane)  but  no 
experimental  measurements  exist.  These  two  bending  modes  are  evidently  very  weak 
and/or  broadened  in  the  argon  matrix  isolation  spectra  as  no  absorptions  positively 
assignable  to  the  dimer  have  yet  been  identified  outside  the  HF  stretching  spectral 
region.  The  tendency  for  the  dimer  H atoms  to  tunnel  and  the  overall  molecular 
f legibility  evidently  affect  the  appearance  of  the  bending  transitions,  most 
likely  broadening  them.  This  type  of  behavior  was  recently  observed  for  certain 
vibrational  modes  of  tropolone(22) , where  sharp  bands  were  observed  upon  deuter- 
ation  of  the  labile  H atom  and  a neon  matrix,  which  is  less  Interactive  than  argon 
was  used.  The  two  lowest  frequency  vibrational  modes  of  the  dimer  are  the 
F**«F  stretch  and  a rocking  of  the  free  HF  moiety.  Both  of  these  frequencies  are 


expected  to  fall  below  the  range  of  the  presently  observed  spectral  region.  They 
are  calculated (21)  to  occur  at  226  and  171  cm  ^ . 

Trimer . No  experimental  evidence  for  the  existence  of  the  H^F^  trimer  is 
present  in  the  literature.  The  various  studies  (vida  infra  for  discussion)  indi- 


cate that  the  dimer,  tetramer,  and  hexamer  have  greater  stability  in  the  vapor 
phase.  However,  trimer  can  exist  in  vapor,  as  suggested  by  the  mass  spec trorae trie 
observation  of  a (nonpolar)  ion  containing  three  F a toms (11). 

The  matrix-isolation  experiments  show  the  presence  of  four  intense  oligomer 
absorptions  in  the  HF  stretching  region.  The  388  1 cm  ^ band  is  due  to  dimer  and 


are  very  closely  associated  as  shown,  for  example,  by  their  parallel  behavior 
in  the  presence  of  SiF^  impurity.  Plots  of  the  three  band  intensities,  scaled 
relative  to  the  normalized  3881  cm  1 dimer  band,  show  that  all  three  increase 
relative  to  the  dimer  band  (but  noL  bv  the  same  amount)  as  the  M/S  ratio  is 
made  more  concentrated.  All  three  bands  are  assigned  to  the  chain  trimer  in 
preference  to  including  cyclic  trimer  or  higher  oligomers.  The  chain  trimer 
is  chosen  because  it  has  three  infrared  active  HF  stretching  modes,  the  inten- 
sities of  the  three  observed  bands  correlate  reasonably  well  together,  and  because 
they  occur  in  the  spectral  region  of  die  (chain)  dimer  HF  stretching  modes.  There 
is  a gap  in  the  appearance  of  the  intense  HF  stretching  modes  in  the  argon  matrix 
spectra  that  strongly  suggests  the  presence  of  chain  dimer  and  trimer  with  absorp- 
tions in  the  range  3702-3920  cm  ^ , and  other  oligomers  with  HF  stretching  absorp- 
tions in  the  range  below  3430  cm  It  is  known  that  cyclic  hexamer  and  (proba- 
bly cyclic)  tetramer  have  their  single  infrared  allowed  mode  in  the  latter 
region(23)  and  that  the  long  chain  (solid)  HF  system  also  absorbs  there(24). 

The  chain  trimer  fits  very  well  into  a double  substitutional  site  in  solid 
argon.  Matrix  warmup  studies  suggest  that  it  is  matrix  stabilized  in  its  trap- 
ping site,  as  the  sharp  trimer  absorption  bands  persist  while  absorptions  due 
to  other  species  disappear.  This  type  of  behavior  has  also  been  observed  for 
LiRFn  oligomers ( 25) . 

Several  weak  absorptions  are  observed  in  the  HF  stretching  region  above 
3700  cm  ^ that  may  be  associated  with  chain  oligomers  longer  than  trimer.  All 
n HF  vibrations  of  the  chain  oligomers  are  infrared  active,  in  contrast  with 
the  single  allowed  mode  of  the  planar,  cyclic  molecules. 

The  chain  trimer  possesses  four  FH*  • • F bending  modes.  Assuming  linear 
hydrogen  bonds,  which  define  a plane,  the  in-phase,  out-of-plane  mode  should 
be  intense  while  the  out-of-phase  mode  should  be  weaker.  The  two  in-plane  modes 


should  be  far  less  affected  by  tunneling  and,  therefore,  nonrigidity  than 


is  the  corresponding  dimer  mode  since  tunneling  involves  motion  by  three  H atoms 
rather  than  two.  In  contrast  with  the  unobserved  dimer  vibrations,  intense 
trimer  bands  attributed  to  FH*  • -F  bending  modes  are  observed  at  446  and  561  cm 
The  remaining  five  chain  trimer  modes  (two  F*  • *F  stretches,  two  "rocking" 
modes,  and  one  torsion)  are  all  expected  to  fall  below  400  era 

Hexamer  and  Tetramer.  The  H^F^  hexamer  is  known  to  exist  abundantly  in  the 
vapor  at  pressures  approaching  saturation  (see  discussion  in  the  next  section), 
but  it  is  most  likely  not  present  in  the  low  pressure  HF  vapor  deposited  in  pre- 
mixed matrix-isolation  samples.  Electron  diffraction  studies (26)  show  that  the 
hexamer  is  cyclic;  however,  it  is  not  positive  that  the  molecule  is  planar. 

Very  likely  it  is  planar  with  large  amplitude  angle  bending  vibrations,  but  the 
possibility  for  boat  or  chair  equilibrium  configurations  also  exists(26). 

A planar  cyclic  hexamer  (symmetry  group  C^)  and  tetramer  (C^)  each  have 
only  four  infrared  active  vibrations.  These  are  an  HF  stretch,  an  in-plane  and 
an  out-of-plane  angle  deformation,  and  a low-frequency  F*  • *F  stretch.  Spectral 
studies  by  Smith(23)  yield  rather  approximate  values  for  the  overlapped,  infrared 
active  HF  stretches  in  the  vapor.  He  gives  derived  band  profiles  with  maxima 
near  3350  and  3500  cm  \ respectively,  for  hexamer  and  tetramt»r.  The  absorption 
is  assigned  to  these  oligomers  on  the  basis  of  its  temperature  and  pressure  depen- 
dence. Smith  reports  additional  absorption  near  1200  and  700-850  cm  ^ for  hexa- 
mer and  near  700  cm  ^ for  tetramer.  The  argon  matrix-isolation  results  have  not 
yet  been  completed  on  these  oligomers  because  the  double  jet  deposition  experi- 
ments have  not  been  perfected.  Nevertheless,  bands  due  to  HF  oligomers  occur 
at  3430  cm  ^ and  3250  cm  ^ which  are  very  likely  due  to  the  tetramer  and  hexamer 


because  of  their  proximity  to  the  estimated  vapor  phase  values.  These  bands  grow 
on  diffusing  dilute  matrix-isolation  samples,  and  occur  only  for  initial  deposits 
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with  concentrated  M/S  ratios. 

The  relative  sparclty  of  bands  in  the  FH-  • *F  deformation  region,  as  well 
as  in  the  HF  stretching  region,  corroborates  the  idea  that  the  larger  oligomers 
of  HF  are  cyclic  and  symmetric  rather  than  chains  and  asymmetric.  The  1200  cm  ^ 
absorption  reported  by  Smith  for  the  vapor  phase  hexamer  was  not  clearly  observed 
in  the  present  experiments;  frequencies  proposed  for  the  single  infrared  active 
angle  deformation  modes  of  the  tetramer  and  hexamer  are  listed  in  the  Tables. 


Hydrogen  Fluoride  - Water  Complexes.  The  FH* 


•OH^  dimer  has  been  studied 


by  microwave  spectroscopy (27) , infrared  spectroscopy (28) , and  ab  initio  quantum 
mechanics (29) . The  structure  is  analogous  to  the  H2F2  dimer,  with  a roughly  linear 


FH*  • *0  linkage  but  it  has  not  been  established  that  it  is  a nonplanar  mole- 
cule. The  published  infrared  spectra(28)  are  rather  incomplete,  with  only  the 
direct  observation  of  a broad  vapor  phase  band  extending  downwards  in  frequency 
from  about  3600  cm  ^ and  another  broad  band  near  650  cm  Hot  bands  are  super- 
posed on  these  absorptions,  and  Thomas (28)  interprets  the  peaks  in  a manner  lead- 
ing to  values  for  four  fundamental  bending  modes.  No  information  exists  concern- 
ing higher  oligomers. 

Argon  matrix-isolation  experiments  using  HF  and  H20  mixtures  yield  absorptions 
clearly  due  to  complexes.  The  available  data  suggest  that  a trimolecular  complex 
containing  two  HF  molecules  and  one  H2O  is  particularly  stable  in  the  argon  matrix. 
A strong  absorption  band  at  3554  cm  ^ grows  into  the  spectrum  and  the  3881  cm  ^ 
h2f2  ^imer  band  concomitantly  disappears  on  concentrating  H20  in  the  matrix.  The 
3554  cm  1 frequency  is  significantly  lower  than  the  3608  cm  1 vapor  phase  value 
given  for  the  FH«  • *0H2  dimer(28),  whereas  the  vapor  phase  and  argon  matrix 
frequency  values  for  the  corresponding  HF-  • -HF  vibration  agree  well.  On  review- 
ing the  argon  matrix  data,  it  appears  that  the  FH*  • -FH*  • *OH0  molecules,  like- 
wise, behave  similarly.  The  dimer  spectra  are  weak  in  relatively  concentrated 
matrices,  while  the  trimer  spectra  show  several  veiy  Intense  bands.  The  trimer 
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molecules  both  fit  double  substitutional  sites  in  argon  and  warmup  studies  sug- 
gest that  both  are  matrix-stabilized.  The  second  HF  stretching  mode  of  the  mixed 
(HF^Oi^  trimer  is  weak  as  it  is  basically  an  antisymmetric  combination  of  the 
two  HF  oscillators.  All  three  "^O"  modes  are  weak  in  its  complexes,  and  other 
research  shows  they  are  not  shifted  much  in  frequency  from  the  monomeric  H?0 
values (9, 30) . 

Four  bands  appear  just  below  650  cm  ^ in  the  FH*  • • F angle  deformation  region. 
These  bending  modes,  which  are  assigned  in  the  Tables,  are  higher  in  frequency 
than  those  of  the  C1H-  • *(^2  dimer (30)  but  close  to  the  vapor  phase  frequencies 
given  for  FH‘  • ’OH^  dimer  by  Thomas(28). 

B.  Composition  of  Hydrogen  Fluoride  Vapor 

It  is  well  known  that  at  modest  temperatures  HF  vapor  is  highly  associated 
at  pressures  extending  considerably  below  its  saturation  vapor  pressure.  In 
contrast  to  most  other  substances,  hydrogen  fluoride  has  a tendency  to  form  hydro- 
gen bonded  vapor  phase  oligomers,  rather  than  to  immediately  condense  to  form  the 
liquid.  In  addition  to  association,  it  is  expected  that  the  strong  attractive 
forces  existing  between  HF  molecules  should  lead  to  significant  nonideality  of 
the  vapor.  An  attempt  is  made  here  to  provide  a unified  model  for  the  HF  vapor 
by  introducing  gas  nonideality  in  the  form  of  a virial  coefficient.  The  nonideal- 
ity contribution  nearly  vanishes  for  low  pressure,  high  temperature  data  sets, 
but  plays  an  important  role  for  data  sets  obtained  at  pressures  near  saturation. 

It  seems  impossible  to  precisely  represent  the  comp  ^x  overall  experimental  phe- 
nomena observed  for  HF  vapor  In  terms  of  very  few  adjustable  parameters;  however, 
a moderately  successful  attempt  is  made  to  reconcile  the  vapor  density  data, 
specific  heat  data,  infrared  spectroscopic  data,  and  third  law  entropy  data  in 
terms  of  one  or  two  adjustable  gas  nonideality  parameters  and  two  or  three  equili- 
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brium  constants. 

The  density  of  HF  vapor  was  analyzed  in  terms  of  a monomer-hexamer  equili- 
brium many  years  ago (31).  Saturated  and  near  saturated  vapor  yields  a linear 
van't  Hoff  plot  over  a 100°C  temperature  range(31,32)  for  the  hexamer ization 
reaction  6 HF  ^ H^F^.  **  therefore  appears  as  a consequence  that  any  model  for 

HF  vapor  composition  must  lead  to  hexamer  as  the  dominant  oligomer  along  the 
saturation  vapor  pressure  line.  Careful  analysis  by  Maclean,  Rossotti  and  Ros- 
sotti(34),  hereafter  MRR,  of  highly  accurate  vapor  density  data  obtained  by  Stroh- 
meier  and  Briegleb(33) , suggests  that  a monomer-dimer-hexamer  equilibrium  fits 
the  data  best.  The  data  covers  pressure  ranges  that  are  considerably  below  the 
saturation  values,  and  the  data  fit  involves  only  low  pressures  of  dimer  and 
hexamer.  The  MRR  analysis  is  based  on  ideal  gas  behavior  and  the  model  fails  to 
reproduce  the  vapor  density  data  as  pressure  increases  and  oligomer  concentrations 
rise  to  high  levels.  Extrapolation  of  their  results  to  the  saturation  vapor  pres- 
sure predicts  association  factors  Z(Z  = apparent  molecular  weight  of  associated 
vapor  divided  by  molecular  weight  of  monomer)  that  are  far  below  those  observed (32) 
for  the  saturated  vapor. 

Infrared  spectral  studies  of  HF  vapor  by  Smith (12, 23)  suggest  that  dimer, 
tetramer,  and  hexamer  are  all  present  in  HF  vapor  in  the  P-T  regime  investigated 
by  Strohmeier  and  Briegleb (33) . The  infrared  results  suggest  that  the  tetramer 
and  hexamer  have  large  absorption  coefficients  and  that  these  molecules  are 
easily  detected  under  P-V-T  conditions  where  their  presence  has  little  or  no 
effect  on  vapor  density  measurements.  The  hexamer,  for  example,  shows  an  intense 
infrared  absorption  at  -77°C  at  a total  HF  pressure  of  only  2-3  torr(23). 

Third  law  entropy  measurements (35)  of  the  saturated  vapor  strongly  reflect 
molecular  association.  On  using  the  association  factor  Z that  is  obtained  from 
vapor  density  measurements,  a monomer-hexamer  model  fails  to  account  for  the 
entropy  by  about  40  J/mol*deg  (the  calculated  vapor  entropy  is  far  below  that 
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of  the  liquid!);  therefore,  the  original  analysis  of  the  entropy  data  used  an 
extended  oligomer  model (35). 

Vapor  phase  heat  capacity  measurements (36)  show  enormous  maxima  that  roughly 
parallel  the  saturation  pressure  curve.  The  maxima  are  unexcelled  among  vapor 
phase  heat  capacities,  and  arise  because  of  the  enthalpies  of  the  various  oli- 
gomerization reactions.  The  Cp  curves  are  dominated  by  a monomer-cyclic  hexa- 
mer  equilibrium,  but  a more  quantitative  fit  is  obtained  upon  introducing  a 
vapor  model  including  a series  of  chain  oligomers (36) . 

In  general,  each  data  set  has  been  analyzed  in  terms  of  a specific  vapor 
model  and  the  various  models  are  not  consistent  with  one  another.  However,  it 
is  apparent  that  the  entropy  and  heat  capacity  data  can  be  approximately  fitted 
by  a monomer-dimer-hexamer  model  if  the  equilibrium  constant  is  lowered  from 
that  calculated  using  the  observed  saturated  vapor  density  data.  Adjustment 
of  K2  and  to  fit  more  extensive  data  is  accomplished  by  introducing  nonideal 
gas  behavior  as  described  below.  Maclean,  Rossotti  and  Rossotti(34)  found  mini- 


mal evidence  for  any  single  additional  oligomer  beyond  dimer  and  hexamer,  although 
(obviously)  an  exact  data  fit  can  be  accomplished  by  adding  several  oligomers 
because  of  the  rather  indeterminate  equilibrium  constants.  Since  tetramer  was 
observed  at  low  pressures  in  the  vapor  phase  infrared  measurements  along  with 
the  hexamer.  the  present  analysis  includes  dimer,  tetramer,  and  hexamer  as  the 
most  important  oligomer  species.  The  analysis  essentially  reproduces  that  of 
MRR  at  the  ideal  gas  level,  but  extension  by  the  introduction  of  a second  virial 
coefficient  leads  to  calculated  values  in  reasonably  close  agreement  with  each 
of  the  four  different  data  sets. 

In  order  to  orient  the  P-T  regimes  of  the  various  data  sets  relative  to 
one  another,  they  are  presented  together  in  Fig.  2.  Rather  than  plotting  iso- 
therms of  Z versus  pressure,  a plot  of  the  excess  volume  versus  pressure,  as 
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origlnally  drawn  by  Vender zee  and  Rodenberg (37) , 19  presented.  The  excess 
volume,  V^,  Is  simply  the  difference  between  the  actual  volume  of  one  gram 
formula  weight  of  associated  HF  vapor  and  the  ideal  gas  volume.  Horizontal  iso- 
therms reflect  the  presence  of  monomer  and  dimer,  while  sharply  dipping  isotherms 
reflect  the  occurrence  of  higher  oligomerizations.  Superposed  on  the  excess  vol- 
ume isotherms  at  appropriate  pressures  are  the  infrared  and  calorimetric  measure- 
ments, as  described  in  the  figure  caption. 

NONIDEAL  VAPOR  MODEL.  Vapor  density  data  are  expressed  using  the  equation 
of  state  for  one  formula  weight  (20.01  g)  of  HF  as 


PV  = RT/Z 


This  equation  is  empirical  and  in  the  previous  work  Z was  attributed  exclusively 
to  molecular  association.  In  the  present  case  It  is  generalized  to 
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where  Z^  is  the  equation  of  state  correction  due  to  association  and  Z^  is  that 
due  to  gas  nonideality  (including  effects  from  all  vapor  phase  species).  Intro- 
ducing a similarly  comprehensive  virial  coefficient  B,  the  equation  of  state  is 


written 


Since  n - 1/ZA> 


PV  * nRT  (1  + nB/V) 


1/Z  - 1/ZAZV  * n (1  4-  nB/V). 
A N 


1/Z  - 1 + nB/V 

n 


1 + zbp/rtza 


gives  the  nonideality  correction.  The  second  virial  coefficient  for  a gaseous 
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mixture  is  written(38) 

B - l B X X (6) 

ij  ] 3 

where  X.  is  the  mole  fraction  of  the  Mi"th  oligomer,  now  represented  as  H F 

i n, 

Calculations  using  simple  model  potential  functions  suggest  that  the  second 
virial  coefficient  for  each  type  of  molecule  will  be  roughly  proportional  to 
its  volume  and,  with  strong  attractive  forces,  to  an  exponential  term: 
exp  (Ae^/RT),  where  Ae^  is  the  depth  of  the  potential  well  resulting  from 
interaction  between  the  n^  and  n^  oligomers.  Because  all  the  . virial  coef- 
ficients cannot  be  realistically  determined  for  the  HF  oligomers,  simplifying 
assumptions  of  the  form 


r \ 

AW 
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were  made.  The  F^  values  (excluding  F^,  always  unity)  were  assumed  to  take  the 

2 

forms  (A)  F^  = s n^,  or  (B)  F^  « s n^  , with  s a scale  factor  given  the  same  value 

for  all  oligomers.  With  F^  = the  assumption  corresponds  to  saying,  for  example, 

that  the  hexamer  is  six  times  as  voluminous  as  the  monomer,  and  that  the  stronger 

(rotationally  averaged)  dipole  and  hydrogen  bonding  attractions  of  the  monomer 

are  approximately  balanced  by  the  stronger  and  more  orientationally  uniform 

attractive  van  der  Waals  forces  of  the  hexamer  (i.e.  Ae^  is  equal  for  the  two 

2 

oligomers).  For  F^  = n^  the  volume  relation  is  still  the  same  but  exp(Ac^^/RT) 

is  six  times  greater  than  exp(Ae, . /RT) . This  implies  that  (Ac,.  - Ae,,)  ^ 

11  oo  11 

4.7  kJ/mol . 

The  computations  tested  assumptions  (A)  and  (B).  First,  the  vapor  density 
data  of  Briegleb  and  Strohmeier (33) , which  is  given  for  five  temperatures  between 
26  and  56°C,  was  fitted  to  obtain  values  for  the  K0  and  K.  or  the  K0,  K,  and  K. 
equilibrium  constants  as  a function  of  the  second  virial  coefficient  B (i.e. 
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B^  and  s) . Rather  than  using  a graphical  method  as  done  by  Maclean,  Rossotti 
and  Rossotti (34) , a simultaneous  least  squares  fit  of  the  vapor  density  data  to 
the  total  pressure  and  to  the  function 

A » l n.  K.  P ni  (8) 

i 1 1 

was  obtained.  The  monomer  pressure  was  obtained  by  numerical  integration  of 
the  Bjerrum  integral,  and  the  results  for  B = 0 agreed  reasonably  well  witli  the 
results  of  MRR,  as  shown  in  Fig.  3.  Unweighted  least  squares  were  used  to  fit 
the  same  data  points  as  used  by  MRR,  and  a 1/Z  weighting  was  used  when  the  Stroh- 
meier-Briegleb  data  was  supplemented  by  data  listed  in  the  paper  of  Vanderzee  and 
Rodenburg (37) . The  latter  authors  include  unpublished  data  in  their  comprehen- 
sive, and  smoothed,  compilation  of  the  thermodynamic  excess  properties  of  hydro- 
gen fluoride  vapor.  It  must  be  noted  that  their  recommended  values  for  the  asso- 
ciation factor  Z along  the  saturation  curve  differ  from  those  of  Jarry  and  Davis (32) 
by  several  percent  at  26°C.  The  1/Z  weighting  generally  favors  the  Strohmeier- 
Briegleb  data  by  a factor  of  2 or  3 over  data  closer  to  the  saturation  curve. 

The  quality  of  the  equilibrium  constants  and  one  or  two  second  virial  coef- 
ficient parameters  were  then  assessed  by  calculating  the  oligomerization  enthal- 
pies, the  third  law  entropies,  and  the  vapor  phase  heat  capacities.  Their  recom- 
mended values  for  these  quantities  are  tabulated  by  Vanderzee  and  Rodenburg ( 37) . 

RESULTS . Since  the  procedure  used  is  semi empirical  it  was  necessary  to  per- 
form many  trial  computations  to  assess  the  effectiveness  of  the  various  vapor 
phase  models.  The  results  of  one  of  these,  which  yields  rather  good  results,  is 
given  in  Figure  3 and  Table  III.  It  is  a monomer-dimer-tetramer-hexamer  model; 
however,  the  results  are  nearly  identical  to  those  with  tetramer  excluded.  In 
general,  including  tetramer  is  ineffectual  for  significantly  improving  the  fit 
to  the  low  pressure  vapor  density  data,  as  previously  concluded  by  MRR,  and  inclu- 
ding tetramer  to  improve  the  fit  to  the  higher  pressure  data  leads  to  an  improved 


-20- 


overall  average  fit,  but  to  a significantly  poorer  fit  to  the  (more  reliable) 
low  pressure  data. 

In  the  present  case,  the  Strohmeier-Briegleb  vapor  density  data  was  very 

closely  fitted  to  obtain  Kj*  and  values.  The  data  were  fitted  using 

unweighted  least-squares , as  mentioned  above,  rather  than  by  the  visual  curve 

fitting  method  of  MRR.  The  two  ideal  gas  fittings  are  compared  in  the  lower 

two  lines  of  Fig.  3.  The  present  AH^  value,  -32.0  kJ/mole,  is  slightly  higher 

than  MRR's  value  of  -27.2  kJ/mole.  The  latter  is  in  the  range  usually  quoted 

for  the  dimerization  energy.  The  AH.  values  are  similar,  the  present  ideal  gas 

b 

value  being  -173  kJ/mole  while  the  MRR  value  is  -169  kJ/mole.  AH,  changes  very 

o 

little  on  introducing  gas  nonideality;  however,  for  given  models  AH£  may  increase 
to  unreasonable  values  when  the  higher  pressure  vapor  density  data  is  fitted. 
Thus,  the  present  nonideal  gas  model  yields  low  pressure  agreement  plus  agreement 
between  the  calculated  and  the  Jarry-Davis  experimental  Z values  for  the  top  set 
of  points  drawn  in  Fig.  3.  Models  giving  very  good  fits  for  the  saturation  Z 
values  have  so  far  given  unsatisfactory  fits  for  the  entropy  and  heat  capacity. 

The  present  illustrative  nonideal  gas  model  is  for  case  (A)  with  s = ] , 

AH^  is  kept  very  near  the  ideal  gas  value  by  adding  :onstant 
increments  of  -0.10  and  -0.15  to  the  ideal  gas  values  of  -log  l^.  The  resulting 
and  values  are  those  yielding  least-squares  fits  to  the  low  pressure  vapor 
density  data  with  the  new  values.  A value  arises  for  each  temperature, 
but  the  five  B^  values  can  be  correlated  (or  generated)  quite  well  with  a two 
parameter  temperature  dependent  function.  The  set  of  equilibrium  constants  were 
then  used  to  calculate  the  results  shown  in  Table  III,  where  it  is  seen  that  the 
third  law  entropy  values  are  about  2 J/mole*deg  too  low  (rather  than  40!),  that 
the  0^  values  have  approximately  the  correct  maximum  value,  but  that  they  peak 
slightly  too  soon.  Allowing  both  B and  s to  be  temperature  dependent  will  improve 
the  fit.  The  present  Z^  fit  is  best  using  a point  intermediate  between  the 
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nonideality  shifts  of  -0.10  and  -0.15  that  were  chosen  for  illustration  and  to 
indicate  sensitivity.  It  is  seen  that  the  pressures  P.  and  P.  are  of  the  same 

4 O 

order  of  magnitude  in  the  P-T  regime  chosen  by  Smith(23)  for  his  infrared  studies, 
so  that  it  is  reasonable  that  both  oligomers  were  observed.  His  oligomerization 
enthalpies  are  near  the  presently  developed  values. 

The  second  virial  coefficients  represented  by  the  values  increase  reason- 
ably as  the  temperature  is  lowered  from  56  to  26°C.  Before  this  investigation  is 
completed  the  value  of  (i.e.  the  monomer  second  virial  coefficient)  will  be 
computed  using  the  program  LINEAR  obtained  from  OCPE(39)  and  the  HF  monomer-mono- 
mer interaction  potential  function  published  by  Jorgensen  and  Cournoyer (40) . 

C.  Ab  Initio  Calculations  for  HnFn  Oligomers 

In  order  to  augment  the  infrared  experiments,  a theoretical  study  of  the 
higher  HnFn  oligomers  has  been  undertaken.  These  studies  are  similar  in  scope 
to  previous  computational  research  on  the  H2F2  dimer (21, 41)  and  make  use  of  the 
Gaussian  70  computer  program  from  QCPE(42).  With  the  exception  of  one  short  note 
that  considers  H3F3  trimer(43),  there  appear  to  be  no  published  calculations  for 
the  higher  HnF n oligomers . 

When  completed,  the  present  study  will  include  geometry  optimization  of  the 
cyclic  and  chain  configurations  for  the  trimer,  tetramer,  pentamer  and  hexamer. 
Since  the  results  are  dependent  on  the  basis  set,  the  optimum  geometry  will  be 
calculated  using  both  minimal  and  split  valence  basis  sets.  The  analysis  includes 
a calculation  of  force  constants  and  vibrational  frequencies.  The  geometry  opti- 
mization utilizes  the  algorithm  of  Payne(44).  At  present,  minimal  basis  set 
(ST0-3G)  calculations  for  the  chain  and  cyclic  trimer  and  tetramer  are  reported. 

In  each  case  the  oligomer  was  considered  as  one  molecule  and  the  total  energy 
was  minimized  with  respect  to  all  of  the  geometrical  parameters. 
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Trimer . Three  parameters  are  necessary  to  describe  the  geometry  of  the 
planar,  cyclic  trimer:  the  F-F  bond  length,  the  F-H  bond  length,  and  the  angle 
between  FH  and  FF  bonds,  called  a.  The  cyclic  trimer  was  first  assumed  to  be 
planar  and  to  have  a axis  in  order  to  optimize  the  bond  lengths  and  a. 

Planarity  was  confirmed  by  calculating  energies  for  rotation  of  the  protons  out 
of  the  plane  into  both  staggered  and  symmetrical  configurations.  The  optimized 
geometrical  parameters,  a figure  giving  coordinate  definitions,  and  the  total 
energy  are  given  in  Table  IV.  Likewise,  the  chain  trimer  was  initially  assumed 
to  be  planar  and  the  five  bond  lengths  optimized  together  and  the  four  in-plane 
angles  optimized  together.  Planarity  was  confirmed  by  energy  calculations  with 
H atoms  moved  out  of  the  plane.  The  energy  was  lowest  for  the  planar  configura- 
tion given  in  Table  IV,  The  barrier  hindering  rotation  of  the  "free"  H atom 
about  the  adjacent  H*  • *F  axis  is  about  3 kJ . 

The  interaction  energies  were  obtained  using  monomer  and  dimer  energies 
already  given  by  Curtiss  and  Pople(21).  They  calculated  a total  energy  of 
-98.57285  hartrees  for  the  monomer  with  an  HF  distance  of  9.56  nm  using  an 
ST0-3G  basis  set.  For  the  dimer,  the  HF  distances  are  9.54  and  9.53  nm  with 
an  energy  of  -197.15447  hartrees.  These  values  yield  an  interaction  energy  of 
-70.5  kJ,  or  -23.5  kJ  per  hydrogen  bond,  when  three  monomers  combine  to  give 
the  cyclic  trimer.  The  interaction  energy  of  a monomer  and  dimer  to  form  cyclic 
trimer  is  -47.4  kJ  or  -23.7  kJ  per  hydrogen  bond  formed.  Corresponding  values 
for  the  chain  trimer  are  -48.5  kJ  and  -25.4  kJ.  At  -24,3  and  -25.4  kJ , the 
ST0-3G  calculation  suggests  slightly  stronger  hydrogen  bonds  for  the  chain 
trimer,  but  the  cyclic  trimer  would  be  more  stable  due  to  the  extra  hydrogen 
bond . 

Tetramer . Like  the  planar  cyclic  trimer,  only  three  geometrical  parameters 
are  required  for  cyclic  tetramer.  This  molecule  was  initially  assumed  to  be  planar 
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and  the  two  bond  lengths  and  angle  a were  optimized.  Then,  planarity  was  con- 
firmed as  the  H atoms  were  rotated  out  of  the  plane  in  both  staggered  and  sym- 
metrical configurations.  Optimized  values  for  the  parameters  and  the  total  energy 
are  listed  in  Table  IV.  The  binding  energy  of  four  monomers  to  give  the  cyclic 
tetramer  is  -185.3  kJ  or  -46.7  kJ  per  hydrogen  bond.  This  value  appears  to  be 
much  too  high,  just  as  the  FF  bond  distance  appears  too  short.  Results  of  the 
calculations  on  the  chain  tetramer  are  included  in  Table  IV,  where  they  may  be 
compared  with  the  chain  trimer. 

The  minimal  ST0-3G  basis  set  calculations  suggest  that  the  cyclic  trimer 
and  tetramer  are  more  stable  than  the  chain  forms.  Experimentally,  the  chain 
trimer  appears  to  be  the  more  stable  in  an  argon  matrix,  perhaps  because  of 
stabilizing  matrix  interactions.  It  is  premature  to  discuss  the  present  results 
in  detail  in  view  of  the  expected  limitations  of  the  minimal  basis  set;  however, 
they  certainly  anticipate  the  stronger  hydrogen  bonding  that  arises  for  the 
cyclic  hexamer  and  long  chains  of  the  solid.  As  pointed  out  above,  there  is  a 
gap  in  the  frequency  range  of  the  HF  stretches  for  short  chains  and  for  the 
cyclic  and  long  chain  9ystem9.  In  addition,  there  is  a pronounced  shortening 
in  the  experimental  FF  bond  lengths  as  seen  on  comparing  the  dimer,  27.9  nm(ll), 
cyclic  hexamer,  25.5  nm(26),  and  solid,  24.9  nm(45).  Of  course,  the  unstrained 
cyclic  systems  have  properties  in  common  with  the  infinite  chain. 

One  objective  in  performing  the  ab  initio  calculations  is  to  assist  in  the 
vibrational  analysis  of  the  observed  spectra  by  providing  relative  ranges  and 
frequency  orderings  for  the  various  oligomers. 

D.  Summary  of  Results 

1.  Infrared  spectra  of  several  HF  and  HF/H^O  mixtures  suspended  in  solid 
argon  were  obtained  in  the  range  400-4000  cm  ^ . 

I 

l—  J 


-24- 


2.  The  data  show  HF  stretching  bands  attributed  to  monomer,  chain  dimer, 
chain  trimer,  cyclic  tetramer  and  cyclic  hexamer.  The  data  show  FH • • *H  angle 
deformation  modes  attributed  to  these  same  oligomers;  however,  not  all  of  the 
anticipated  bending  frequencies  were  located  (in  particular,  the  dimer  bending 
modes  are  not  listed). 

3.  There  seems  to  be  parallel  behavior  for  the  dimers  FH*  • *FH  and 

FH*  • *FH*  • The  trimer  molecules  appear  to  be  stabilized  in  the  argon 

matrix,  and  to  possess  some  very  intense,  sharp  absorptions. 

4.  The  dominant  oligomers  in  HF  vapor  appear  to  be  the  dimer  and  the  hexa- 
mer. The  latter  molecule  is  very  important  as  saturation  vapor  pressures  are 
approached.  The  next  most  important  oligomer  is  the  tetramer,  although  it  plays 
a very  minor  role.  These  conclusions  are  suggested  by  a literature  review  and 
by  application  of  a nonideal  gas  model  to  vapor  phase  density,  third  law  entropy, 
heat  capacity,  and  infrared  spectral  data. 

5.  There  appears  to  be  stabilization  of  the  FH*  * *F  hydrogen  bonds  on 
catenation  or  on  unstrained  cyclization  of  HF  systems.  Minimal  basis  ST0-3G 
calculations  on  various  oligomers  suggest  this  trend,  which  is  apparent  in  the 
experimental  data. 


E.  References 

1.  M.  T.  Bowers,  G.  I.  Kerley  and  W.  H.  Flygare,  J.  Chem.  Phys . , 4_5,  3399  (1966). 

2.  D.  W.  Robinson  and  W.  G.  Von  Holle,  J.  Chem.  Phys.,  44,  410  (1966). 

3.  M.  G.  Mason,  W.  G.  Von  Holle  and  D.  W.  Robinson,  J.  Chem.  Phys.,  _54,  3491  (1971). 

4.  R.  L.  Redington  and  D.  E.  Milligan,  J.  Chem.  Phys.,  39.*  1276  (1963). 

5.  L.  Abouaf-Marguin,  M.  E.  Jaycox  and  D.  E.  Milligan,  J.  Molec.  Spectrosc.,  67 , 

34  (1977). 

6.  L.  F.  Keyser  and  G.  W.  Robinson,  J.  Chem.  Phys.,  44.,  3225  (1966). 


-25- 


7.  A.  Cabana,  G.  B.  Savitsky  and  D.  F.  Homig,  J.  Chem.  Phys.,  3j),  2942  (1963). 

8.  L.-C.  Brunei,  J.-C.  Bureau  and  M.  Peyron,  Chem.  Phys.,  28^,  387  1978). 

9.  G.  P.  Ayers  and  A.  D.  E.  Pullin,  Spectrochim.  Acta,  32A,  1689  (1976); 

L.  Fredin,  B.  Nelander  and  G.  Ribbegard,  J.  Chem.  Phys.,  6*6,  4065  (1977), 

10.  Precision  Aperture,  Fort  Wayne,  Indiana 

11.  T.  R.  Dyke,  B.  J.  Howard  and  W.  Klemperer,  J.  Chem.  Phys.,  _56,  2442  (1972). 

12.  D.  F.  Smith,  J.  Molec.  Spectrosc.,  _3,  473  (1959). 

13.  W.  F.  Herget,  N.  M.  Gailar,  R.  J.  Lovell  and  A.  H.  Nielsen,  J.  Opt.  Soc. 
Amer.,  50,  1264  (1960). 

14.  J.  L.  Himes  and  T.  A.  Wiggins,  J.  Molec.  Spectrosc.,  AC) , 418  (1974). 

15.  T.  Ueda  and  T.  Shimanouchi,  J.  Molec.  Spectrosc.,  _28,  350  (1968). 

16.  F.  A.  Miller  and  W.  G.  Fatelely,  Spectrochim.  Acta,  20,  253  (1964). 

17.  W.  J.  Lafferty,  A.  G.  Maki  and  E.  K,  Plyler,  J.  Chem.  Phys.,  _40,  224  (1964). 

18.  P.  R.  Bunker  and  B.  M.  Landsberg,  J.  Molec.  Spectrosc.,  67,  374  (1977). 

19.  R.  L.  Redington  and  K.  C.  Lin,  J.  Chem.  Phys.,  _54,  4111  (1971). 

20.  Y.  Marechal  and  A.  Witkowski,  J.  Chem.  Phys.,  48,  3697  (1968). 

21.  L.  A.  Curtiss  and  J.  A.  Pople,  J.  Molec.  Spectrosc.,  6^,  1 (1976). 

22.  R.  L.  Redington  and  T.  E.  Redington,  J.  Molec.  Spectrosc.,  in  press. 

23.  D.  F.  Smith,  J.  Chem.  Phys.,  28*  1040  (1958). 

24.  J,  S.  Kittelberger  and  D.  F.  Hornig,  J.  Chem.  Phys.,  46^,  2099  (1967). 

25.  R.  L.  Redington,  J.  Chem.  Phys.,  4j4,  1238  (1966). 

26.  J.  Janzen  and  L.  S.  Bartell,  J.  Chem.  Phys.,  50,  3611  (1969). 

27.  J.  W.  Bevan,  A.  C.  Legon,  D.  J.  Millen  and  S.  C.  Rogers,  J.  Chem.  Soc.  Chem. 
Comm.,  341  (1975). 

28.  R.  K.  Thomas,  Proc.  Roy.  Soc.  Lond.  A.,  344 , 579  (1975). 

29.  D.  G.  Lister  and  P.  Palmieri,  J.  Molec.  Structure,  3^,  295  (1977). 

30.  B.  S.  Ault  and  G.  C.  Pimentel,  J.  Phys.  Chem.,  77,  57  (1973). 


-26- 


31.  R.  W.  Long,  J.  H.  Hildebrand  and  W.  E.  Morrell,  J.  Amer.  Chem.  Soc.,  63, 

182  (1943). 

32.  R.  L.  Jarry  and  W.  Davis,  Jr.,  J.  Amer.  Chem.  Soc.,  5J_,  600  (1953). 

33.  W.  Strohmeier  and  G.  Briegleb,  Z.  Elektrochem. , 51_>  662  (1953). 

34.  J.  N.  Maclean,  F.  J.  C.  Rossotti  and  H.  S.  Rossotti,  J.  Inorg.  Nucl.  Chem., 

24,  1549  (1962). 

35.  J.  H.  Hu,  D.  White  and  H.  L.  Johnston,  J.  Amer.  Chem.  Soc.,  7_5,  1232  (1953). 

36.  E.  U.  Franck  and  F.  Meyer,  Z.  Elektrochem.,  63-,  571  (1959). 

37.  C.  E.  Vanderzee  and  W.  W.  Rodenburg,  J.  Chem.  Thermo dynamics,  2,  461  (1970). 

38.  J.  0.  Hirschf elder , C.  F.  Curtiss  and  R.  B.  Bird,  "Molecular  Theory  of  Gases 

and  Liquids",  John  Wiley  & Sons,  New  York,  N.  Y.  (1954). 

39.  S.  Murad,  Q CPE  10,  357  (1978). 

40.  W.  L.  Jorgensen  and  M.  E.  Cournoyer,  J.  Amer.  Chem.  Soc.,  100,  4942  (1978). 

41.  A.  Karpfen  and  P.  Schuster,  Chem.  Phys.  Lett.,  44_,  459  (1976). 

42.  W.  J.  Hehre,  W.  A.  Lathan,  R.  Ditchfield  and  J.  A.  Pople,  QCPF  11_,  236  (1973). 

43.  M.  Kertesz,  J.  Roller  and  A.  Azman,  J.  Molec.  Struct.,  3J5,  336  (1977). 

44.  P.  W.  Payne,  J.  Chem.  Phys.,  65,  1920  (1976). 

45.  M.  Atoji  and  W.  N.  Lipscomb,  Acta  Cryst.,  7>  173  (1954). 


i 


-27- 
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Each  of  the  three  subsections  discussed  in  Section  II  of  this  report  is 
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titles,  authors  and  journals  are: 

1.  "Infrared  Matrix-Isolation  Studies  of  HnFn  Oligomers"  by  R.  L.  Reding ton 
and  D.  F.  Hamill;  J.  Chem.  Phys.,  J.  Phys.  Chem.  or  J.  Molec.  Spectrosc. 

2.  "Nonideal  Gas  Model  for  HF  Vapor"  by  R.  L.  Redington;  J.  Phys.  Chem.  or 
J.  Am.  Chem.  Soc. 

3.  "Gaussian  Basis  Studies  of  HnFn  Oligomers"  by  D.  F.  Hamill  and  R.  L. 
Redington;  J.  Phys.  Chem.  or  Int.  J.  Quant.  Chem. 

IV.  PERSONNEL 

1.  Principal  Investigator:  Richard  L.  Redington,  Professor  of  Chemistry,  Texas 
Tech  University. 

2.  Research  Assistant:  Delphia  F.  Hamill.  This  research  will  form  a portion 
of  her  Ph.D.  thesis,  to  be  completed  in  1980. 

V.  MEETINGS 

1.  "Studies  of  HnFn  Oligomeric  Species"  by  R.  L.  Redington;  Air  Force  Contractors 
Meeting,  U.  S.  Air  Force  Academy,  Colorado  Springs,  October,  1979. 

2.  "Infrared  Matirx-Isolation  Spectra  of  HnFn  Oligomers"  by  D.  F.  Hamill  and 
R.  L.  Redington;  Southwestern  Regional  Meeting  of  the  American  Chemical 
Society,  Austin,  Texas,  December,  1979. 
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Table  I 

Infrared  Spectrum  of  Hydrogen  Fluoride 
Isolated  In  an  Argon  Matrix 


Ref.  3 

Ref.  1 

Present 

work 

As  s ignme  n t / e omme  n t 

-1 

3997 

3997  cm 

monomer,  R(l) 

3962 

3963.0 

3962 

monomer,  R(0) 

3953 

3954.3 

3953 

monomer 

3934 

3920 

3919.4 

3920 

chain  t rimer 

3915.6 

3904 

3897.6 

3898 

3878 

3881.4 

3881 

dimer 

3875 

3876 

monomer,  P(l) 

3867.9 

3826 

3826 

chain  trimer 

3818 

satellite;  decrease  on  annealing 

3785 

3787 

dimer 

3777.2 

3774 

water  monomer 

3767 

3757 

3757 

water  monomer 

3755.2 

3744.0 

3444 

water  dimer 

3737.7 

(3737) 

(3726) 

3724.4 

3718 

HF-water;  trimer?,  OH  stretch 

3712 

water  monomer 

3711.0 

3702.4 

3702 

chain  trimer 

3697.6 

3694 

satellite;  inc . , then  dec.  on  annea 

3675 

3669.0 

3669 

HF-water;  trimer?,  HF  stretch 

(3601) 


3637.9  cm 
3656 
3606 
3596.4 
3573.3  3573 

3556  3554 

3530  3533 

3517 
3483 

3428 

3414 

3403 

3374 

3250-60 

753 

746 

723 

716 

698 

683 

636 

630 

615 

609 

579 

561 

499 

464 

446 


water  dimer 

HF  - water  trimer 


cyclic  tetramer 


cyclic  hexaraer 

cyclic  hexamer? 
satellite 
cyclic  tetramer? 
satellite 

increase  on  annealing 
increase  on  annealing 
HF-water,  trimer? 
HF-water,  satellite 
HF -water,  trimer? 
HF-water,  satellite 
increase  on  annealing 
chain  trimer 

increase  on  annealing 
chain  trimer 
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Table  III,  Part  A 

* 

Thermodynamic  Properties  Calculated  for  HF  Vapor  Models 


Oligomerization  Enthalpies,  kJ/mole 


Model 

AH2 

AH. 

4 

AH . 

o 

a) 

-32.1 

-88.0 

-172.5 

b) 

-32.0 

-92.8 

-172.1 

Entropy  of  Vapor 

, J /mole* deg 

Model 

56°C 

44°C 

38°C 

32°C 

26°C 

a) 

78.1 

68.2 

63.4 

59.8 

55.9 

b) 

78.4 

68.4 

63.8 

60.1 

56.2 

observed 

78.7 

67.9 

65.4 

63.1 

58.5 

Nonideality  factor,  Zn, 

of  Saturated  Vapor 

Model 

56°C 

44°C 

38°  C 

32°C 

26°C 

(a) 

1.14 

1.20 

1.18 

1.24 

1.30 

observed 

1.35 

1.28 

1.25 

1.26 

1.27 

(b) 

1.21 

1.30 

1.27 

1.39 

1.49 

observed 

1.35 

1.28 

1.26 

1.26 

1.28 

Monomer  Second  Vlrial  Coefficients, 

dm  ^ /mole 

Model 

56°C 

44°C 

38°C 

32°C 

26°C 

(a) 

-0.41 

-0.73 

-0.78 

-1.13 

-1.53 

(b) 

-0.58 

-1.02 

-1.08 

-1.61 

-2.16 

Models  (a)  and  (b)  described  In  text  and  also  shown  in  van't  Hoff  plots  of 
Fig.  3. 
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Table  IV 


Optimized  Geometries  Using  STO-3G 

Basis  Set 

Coordinate 

H3F3  ring 

H4F4  ring 

H3F3  chain 

H.F,  chain 
4 4 

FH,  rx 

9.59  nm 

9.38 

9.55 

9.56 

r2 

9.54 

9.55 

r3 

9.53 

9.57 

r4 

9.52 

FF,  R2 

25 . 8 nm 

23.1 

25.0 
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Figure  2«  Compilation  of  data  showing  nonideality  of  HF  vapor.  The 
dominant  curves  show  the  excess  volume,  V^,  of  HF  vapor  as  a function 
of  pressure.  "Max  C " shows  the  position  of  maxima  in  the  heat  capa- 
city curve.  "MDH"  ^ shows  the  extent  of  highly  accurate  fitting  of 
vapor  density  data  by  a monomer-dimer-hexamer  model,  "MD"  gives  the 
extent  of  the  fitting  by  a monomer- dimer  model.  The  " 0"  dots  indicate 
the  pressure  ranges  of  the  accurate  vapor  density  data;  the  "x"  points 
and  lines  indicate  the  pressure-temperature  regimes  of  the  infrared 
measurements . 
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MRR:  van't  Hoff  plots  for 
the  dimerization  reaction  in 
HF  vapor  phase  models  obtained 
by  fitting  vapor  density  data. 


IDH  and  IDTH:  ideal  gas  results 
for  the  present  least-squares  fitting 
of  the  same  data  with  and  without  tetramer. 


ZO:  nonideal  gas  model  that  fits  both  the 

low  pressure  curve  plus  the  saturated  vapor  density 
density  point. 


(A)  and  (B):  the  IDTH  points  raised  by  -0.10  and  - 0.15 
to  introduce  gas  nonideality  but  maintain  the  ideal  gas 


AH ^ value. 
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Figure  3.  van1 t Hoff  plots  for  the  dimerization  reaction  in  HF 
vapor  phase  models  obtained  by  fitting  vapor  density  data. 
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